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ABSTRACT 

We present an echelle spectrum of the Swift GRB 050730, obtained four hours after the burst using 
the MIKE spectrograph on the Magellan Clay Telescope when the afterglow was at R — 17.7. The 
spectrum reveals a forest of absorption features superimposed on a simple power-law shaped continuum, 
best described as f v {\) oc A" with a = 1.88 ± 0.01 over A = 7000 - 9000 A. We identify the GRB host at 
zgkb = 3.96855 based on the hydrogen Lyman absorption series, narrow absorption lines due to heavy 
ions such as OI, CII, Sill, S II, Nill, Fell, CIV, Si IV, and NV, and fine structure transitions such as 
OF, OF*, Si IF, CIF, and Fell*. Together these transitions allow us to study the the properties of 
the interstellar medium (ISM) in the GRB host. The principal results are as follows. (1) We estimate a 
neutral hydrogen column density of log N(H I) = 22. 15 ±0.05 in the host. (2) The associated metal lines 
exhibit multiple components over a velocity range of ~ 80 km s~ x , with > 90% of the neutral gas confined 
in 20 km s _1 . (3) Comparisons between different ionic transitions show that the host has little/no dust 
depletion and has 1/100 solar metallicity. (4) The absorbing gas has much higher density than that 
of intervening damped Lya absorption (DLA) systems. In addition, we report the identification of an 
intervening DLA system at zdla = 3.56439 with log N(H I) = 20.3 ± 0.1 and < 5% solar metallicity, 
a Lyman limit system at zlls = 3.02209 with log N(H I) = 19.9 ± 0.1, a strong Mgll absorber at 
ZMgii = 2.25313, and a pair of Mgll absorbers at ZMgii = 1-7731, 57 km s -1 apart. We demonstrate 
that rapid echelle spectroscopy of GRB afterglows helps to reveal a wealth of information in the ISM 
and the intergalactic medium along the sightline which, when followed up with late-time, deep imaging, 
will allow us to uncover a sample of distant galaxies with known ISM properties to constrain galaxy 
formation models. 

Subject headings: gamma rays: bursts — ISM: abundances — ISM: kinematics — intergalactic medium 



1. INTRODUCTION 

Various surveys designed to detect emission at optical, 
' near-infrared, and sub-mm wavelengths have yielded large 
, samples of galaxies at redshift z > 2 (e.g. Steidel et al. 
1999; Blain et al. 2002), but whether these galaxies are 
representative of the galaxy population at high redshifts 
and how they are related to the local population is not 
clear. Damped Lya absorption (DLA) systems probe high- 
redshift gaseous clouds of neutral hydrogen column den- 
sity N(H I) consistent with what is observed in the disks 
of nearby luminous galaxies (e.g. Wolfe et al. 2005). They 
are selected uniformly with V(H I) > 2 x 10 20 cm~ 2 and 
represent a unique sample of distant galaxies with known 
interstellar medium (ISM) properties from absorption-line 
studies (e.g. Pettini et al. 1999; Prochaska et al. 2003). 
Indentifying the stellar counterpart of the DLAs has, how- 
ever, been challenging because of the glare of background 
quasars (Colbert & Malkan 2002; Le Brun et al.1997; Rao 
et al. 2003; Chen & Lanzetta 2003). 

Long-duration gamma-ray bursts (GRBs) are believed 
to originate in the death of massive stars (e.g. Woosley 
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1993; Paczyhski 1998; Bloom et al. 2002; Stanek et al. 
2003). Some bursts are followed by optical afterglows (e.g. 
Akerlof et al. 1999) that can briefly exceed the absolute 
brightness of any known quasar by orders of magnitude 
and serve as bright background sources for probing in- 
tervening gas along the line of sight. Because of their 
transient nature, however, optical afterglows do not inter- 
fere with follow-up studies of absorbing galaxies close to 
the sightlines. Early-time, high-resolution spectroscopy of 
GRB afterglows together with deep, late-time imaging of 
galaxies along the sightlines offers a novel means to un- 
cover a sample of high-redshift galaxies based on their ab- 
sorption properties. In the context of heirarchical struc- 
ture formation, GRB progenitors can form before massive 
black holes, and therefore may be used to probe early uni- 
verse, well into the age of reionization. 

Prompt localization of GRB afterglows is critical in or- 
der to take advantage of their brief but extreme brightness 
for acquiring echelle spectroscopy. It has been difficult in 
the past to carry out rapid spectroscopy for the optical 
transient (OT) because of time-consuming processes to lo- 
calize the bursts. Despite an extensive effort, only a small 
number of GRBs have been spectroscopically identified at 
z > 2 with low-to-moderate resolution spectra available 1 
and two with echelle data available (Fiore et al. 2005). 
Together these data show that all GRB host galaxies have 
abundant neutral gas, and some have the largest N(H I) 
among all DLA systems (Jensen et al. 2001; M0ller et al. 

1 see http://www.mpe.mpg.de/~jcg/grbgen.html for a complete 
list. 
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2002; Castro et al. 2003; Jakobsson et al. 2004; Vreeswijk 
et al. 2004). 

The current generation gamma-ray satellite, Swift, is de- 
signed to provide nearly-instant localization of new GRBs. 
Over the past two years, our group has been pursuing well- 
localized GRB afterglows with moderate-to-high resolu- 
tion spectroscopy (e.g. Prochaska et al. 2004; Bloom et al. 
2005). The primary goal of our project is to collect a sta- 
tistically significant sample of galaxies along the sightlines 
toward high-rcdshift afterglows. This galaxy sample in- 
cludes both GRB host galaxies and those foreground DLA 
galaxies that are close to the sightlines. While interven- 
ing DLA systems arise preferentially in the outskirsts of 
distant galaxies (because of a gas cross-section selection 
effect), the GRB host sample offers a unique opportunity 
to study ISM physics more immediate to vigorous star- 
forming regions in high-redshift galaxies (that have rela- 
tively small cross-section). 

Here we present the first echelle spectrum of a Swift 
GRB, obtained four hours after the burst. The spectrum, 
which spans a wavelength range from 3300 A through 9400 
A, exhibits abundant absorption features superimposed on 
a simple power-law shaped continuum. We identify the 
GRB at zqrb = 3.969 based on a strong damped ab- 
sorption trough centered at 6040 A and a suite of asso- 
ciated metal absorption features. In addition, we iden- 
tify an intervening DLA system at zdla = 3.564 and a 
number of strong Mgll absorbers along the sightline. We 
demonstrate that rapid echelle spectroscopy is plausible for 
well-localized afterglows and that the observations provide 
new insight into the nature of GRB host environments, as 
well as the physical properties of the intergalactic medium 
(IGM) along the sightlines. 

2. OBSERVATIONS AND DATA REDUCTION 

We observed the OT of GRB050730 that was first re- 
ported by Hollan et al. (2005) and later confirmed by Sota 
et al. (2005) at RA(J2000) = 14 h 08 m 17M4 and Dec(J2000)= 
-03 d 46'17"8, using the MIKE echelle spectrograph (Bern- 
stein et al. 2003) on the 6.5 m Magellan Clay telescope 
at Las Campanas Observatory. The spectrograph con- 
tains a blue camera and a red camera, allowing a full 
wavelength coverage from near-ultraviolet 3300 A through 
near-infrared 9400 A. The observations were carried out 
in a sequence of three exposures of duration 1800 s each, 
starting at UT 00:00 on 2005 July 31 (four hours after 
the initial burst) when the OT had R « 17.7 (Holman 
et al. 2005). The mean seeing condition over the period 
of integration was 0.7". We used a 0.7" slit and 2x2 
binning during readout, yielding a spectral resolution of 
FWHM « 10 km s" 1 at wavelength A = 4500 A and w 12 
km s _1 at A = 8000 A. The data were processed and re- 
duced using the MIKE data reduction software developed 
by Buries, Prochaska, & Bernstein 2 . Wavelengths were 
calibrated to a ThAr frame obtained immediately after 
each exposure and subsequently corrected to vacuum and 
heliocentric wavelengths. Flux calibration was performed 
using a sensitivity funtion derived from observations of the 
flux standard NGC7293. 

The final stacked spectrum of the afterglow is presented 
in Figure 1. The signal-to-noise of the echelle data is 

2 see http://web.mit.edu/~burles/www/MIKE/mikc_cookbook.html. 



S/N = 12 at 4800 A and S/N = 13 at 8000 A per res- 
olution element. We have performed a \ 2 fitting routine 
over the spectral region between 7000 A and 9000 A, corre- 
sponding to a wavelength range between 1400 A and 1800 
A in the rest frame of the host) and obtained a best-fit 
power-law model f v {\) oc A" with a — 1.88 ± 0.01, which 
is presented as the dotted curve in Figure 1. The best- 
fit power-law index is significantly steeper than what has 
been measured for the OT of GRB020813 (a = 1) at rest- 
frame optical wavelengths (Barth et al. 2003). 

3. ANALYSIS 

The echelle spectrum of the OT of GRB050730 exhibits 
a large number of absorption features due to the H I Lya 
transition and various heavy ions in the ISM of the GRB 
host as well as in intervening gaseous clouds at z < zgrb- 
Most notably is the strong damping trough centered at 
roughly 6040 A, blueward of which we observe a forest of 
absorption features that are not present on the red side of 
the damped absorption profile. We identify the damped 
absorber as the DLA originating in the host galaxy of the 
GRB at zgrb = 3.969 and the forest lines as the Lya forest 
at z < zqrb- In addition, we also identify an intervening 
DLA at zdla = 3.564 and a number of strong metal-line 
absorbers. Here we summarize the physical properties of 
these strong absorbers. 

3.1. The GRB Host Environment at z G rb = 3.969 

The redshift measured for the GRB host using the DLA 
feature is confirmed by associated metal-line transitions. 
We measure a more precise redshift of the GRB at zgrb = 
3.96855 ± 0.00005 using narrower metal absorption lines. 
At this redshift, we measure A^(H I) of the GRB host by 
fitting Voigt profiles to the observed Lya and Ly/3 transi- 
tions simultaneously using the VPFIT software package 3 . 
We obtain log N(U I) = 22.15 ± 0.05 for the host, with 
the error estimated from varying the continuum level of 
the absorption line profiles. This is the highest N(H I) 
observed in DLA systems, including those arising in GRB 
hosts (c.f. Vreeswijk et al. 2004). 

In addition to absorption features due to neutral hydro- 
gen, the spectrum also exhibits a suite of metal-absorption 
lines at the GRB host redshift (Figure 2), from neutral 
species such as O I, to low-ionization transitions such as 
C II, Si II, S II, Ni II, and Fe II, and to high-ionization tran- 
sitions such as C IV, Si IV, and N V. We also identify strong 
fine structure lines such as OF, 01**, Sill*, CII*, and 
Fell*. With the exception of CII* (Wolfe, Prochaska, & 
Gawiser 2003), none of these transitions has been detected 
in intervening DLA systems. In particular, Fell* absorp- 
tion features are only observed locally in Eta Carinae (Gull 
et al. 2005). These saturated fine structure transitions in- 
dicate an extreme ISM environment with high gas density 
that is rarely observed in intervening DLA systems. Fur- 
thermore, the profiles of well resolved lines (e.g. S II 1250) 
show that > 90% of the neutral gas is confined to a veloc- 
ity width 20 km s _1 , which is considerably smaller than 
the median value of intervening DL A systems and implies 
a quiescent environment. But the profiles of saturated 
lines do exhibit absorption extending to « 80 km s -1 (e.g. 
Sill 1526). Finally, the asymmetry of these line profiles 

3 see http://www.ast.cam.ac.uk/~rfc/vpfit.html. 
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Fig. 1.— Stacked echelle spectrum of the OT of GRB050730. The signal-to-noise of the echelle data is S/N = 12 at 4800 A and S/N = 13 
at 8000 A per resolution element. Blueward of the Lya absorption trough at 6040 A from the GRB host environment is a forest of foreground 
Lya absorbers, including a DLA system at 5550 A (zdla = 3.564). The Ly/3 absorption feature of the GRB host is apparent at 5090 A. The 
absence of flux at wavelengths below 4530 A indicates little/no high energy photons beyond the Lyman limit transition escape the host of the 
GRB. The continuum at A = 7000 - 9000 A is best-described by a power-law model f v (X) oc A 1 - 88 as shown in the solid *red) curve. Note 
that the power-law fit from 7000A to 9000A somewhat underpredicts the apparent continuum flux from 4600 A to 5000 A (dotted curve). 



is suggestive of an organized velocity field, e.g. rotation or 
outflow. 

We measure the column density of each transition using 
both the apparent optical depth method (Savage & Sem- 
bach 1991) and the VPFIT software package. Comparing 
column density ratios between different transitions allow 
us to constrain the physical properties of the ISM in the 
GRB host. The principal results are as follows. (1) We 
measure [S/H] = —2.0 ± 0.1. Because S is non-refractory, 
its gas-phase abundance gives a direct measurement of the 
gas metallicity. (2) We find [S/Fc] = +0.3, consistent 
with the gas phase [a/Fe] measurements of low- metallicity 
DLA systems. Even if we adopt an intrinsic solar abun- 
dance pattern, the dust-to-gas ratio in the host ISM is very 
low (c.f. Savaglio et al. 2003). (3) We measure [N/S] = 
— 1.0 ±0.2, again consistent with low-metallicity DLA sys- 
tems (Prochaska et al. 2002). (4) We detect no molecular 
lines in the Lya forest, suggesting a warm gas phase. (5) 
Based on the observed ratio of iV(FeII*)/AT(FeII), we in- 
fer a number density «h > 10 3 cm~ 3 for a temperature 
T < 30000 K. This constrains the size of the host DLA 
cloud to be Zdla < 4.6 pc. 



In summary, aside from the large gas density nn in- 
ferred from saturated fine structure lines, the ISM of the 
GRB host has very similar characteristics to known DLA 
systems at z ~ 4. 

3.2. Intervening Absorbers 

We have identified a number of strong absorbers along 
the sightline toward GRB050730. We briefly summarize 
each system below: 

DLA system at zdla = 3.56439: We identify Lya and 
Ly/3 transition for this absorber and measure log N(H I) = 
20.3 ± 0.1. In addition, we also find associated metal- line 
absorbers due to Sill, A1II, Fell, Si IV, and CIV. The 
weak line strengths of these metal absorption lines suggest 
a low metallicity in the neutral gas, [Si/H] < —1.3. 

Lyman Limit system at zlls = 3.02209: We measure 
N(R I) using the Lya feature and find log N(H I) = 19.9± 
0.1 for this Lyman limit system. In addition, we identify 
Sill, A1II and Fell, but we do not detect CIV absorption 
at the absorber redshift. Including no ionization correc- 
tion, we find [Si/H] = -1.5 ± 0.2. 

Mg II system at ZMgii = 2.25313: We identify saturated 
Mgll doublet at this redshift. In addition, we find strong 
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Fig. 2. — Absorption profiles of different ionic transitions found at the redshift of the GRB host. The l-cr error spectrum is presented as 
the thin (cyan) curve in each panel. The zero relative velocity corresponds to redshift z = 3.96855. 



MgI2852 and Fell absorption features for this absorber. 

Double Mgll systems at ZMgii = 1-7731: We identify a 
pair of Mgll absorbers at Av = 57 km s" 1 apart, for which 
we also find MgI2852 and the Fell absorption series. 

4. DISCUSSION AND CONCLUSIONS 

GRB afterglows clearly provide a novel alternative to 
quasars for probing the ISM and IGM in the young uni- 
verse. In the case of GRB050730, the afterglow reached an 
initial brightness of R = 15.5 (Klotz et al. 2005) and faded 
to R = 17.5 four hours later when the echelle spectroscopy 
was commenced. At z = 3.969, R = 15.5 corresponds to 
an absolute magnitude Mi 450 = —29.3 + 51og/i at rest- 
frame 1450 A 4 . This intrinsic luminosity, within a minute 
after the burst, is comparable to the most luminous quasar 
known (c.f. HS 1700+6416 at z = 2.73; Schneider et al. 
2003). Even at R = 17.5, the corresponding intrinsic lu- 
minosity easily competes with the brightest quasars known 
at this redshift (see e.g. Fan et al. 2001). While absorption 
line systems uncovered toward the sightlincs of optically 
selected quasars are likely to miss dusty absorbers, the 
extreme brightness of early GRB afterlows offers an unbi- 
ased view of the IGM at all epochs. In addition, the simple 
power-law shaped continuum of a GRB afterglow allows a 
more precise and accurate measurement of IGM opacity. 

The sightline toward GRB050730 is particularly inter- 
esting with a strong DLA feature arising in the GRB host 
galaxy and an intervening DLA system at lower redshift. 

4 Adopting a A cosmology, fijn = 0.3 and Oa = 0.7, with a di- 
mensionless Hubble constant h = i^o/(100 km s -1 Mpc ) 



The GRB host DLA system (zdla = z grb), which are 
presumably selected by vigorous star formation and there- 
fore probe deep into the center regions of distant galaxies, 
presents a nice contrast to the intervening DLA system 
(zdla < ^grb), which arises preferentially at large galac- 
tocentric radii due to a larger cross-section of the outskirts 
than the inner regions. Our study shows that aside from 
having a much higher neutral gas density, the GRB host 
DLA has very similar characteristics to known z ~ 4 DLA 
systems, such as low dust content, low metallicity, and 
a-element enhanced chemical composition. 

We have also shown that this sightline runs through a 
number of strong intervening absorbers, including a Ly- 
man limit system and two strong Mgll absorbers. De- 
tailed analyses of various ionic transitions associated with 
these absorbers allows us to study physical properties of 
ISM and IGM over a wide redshift range from z = 1.7 
through z = 3.9. Follow-up deep imaging and low-resolution 
spectroscopy of faint galaxies along the line of sight will 
allow a direct comparison between the physical properties 
of the cold ISM (such as metallicity, kinematics, and dust 
content as derived from absorption line studies) and stel- 
lar properties (such as luminosity, morphology, and star 
formation rate as extracted from absorbing-galaxy anal- 
yses). A large sample of GRB host DLA systems offers 
a unique opportunity to dissect the ISM properties that 
are directly connected to GRBs, while a statistical sample 
of intervening DLA systems identified toward GRB sight- 
lines will provide important insights toward understanding 
the nature of high-redshift DLA systems and offer a sim- 
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pie test to discriminate between different galaxy formation 
scenarios (e.g. Haehnelt, Steinmetz, & Rauch 2000). 
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